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I present a method for calculating heavy quark production which extrapolates 
smoothly from the well understood limits of the fixed-flavour number scheme for 
~ to the zero-mass scheme at Q^/m?jj — > oo. For all > m|j the 

evolution of the heavy quark distribution is precisely as in the massless MS scheme. 
The method is simple to implement and compares well with data. 



1 Introduction 

Until recently the treatment of heavy flavours in global analyses of structure 
functions has been rather naive. Usually the charm and bottom quarks were 
regarded as infinitely massive below = rnj^ , and massless above this thresh- 
old. This procedure, known as the zero mass variable flavour number scheme 
(ZM-VFNS), gives the correct results for Q^/m^ —> oo, but has errors of 
0{m'^jj /Q'^) and is very inaccurate near threshold. Alternatively, all heavy 
flavour flnal states may be thought of as being produced by the hard scatter 
between the electroweak boson and a light parton, i.e. the number of active 
flavours, is 3 and the heavy flavour cross-section is generated (mainly) by 
photon-gluon fusion. This is called the fixed flavour number scheme (FFNS), 
and incorporates the correct threshold behaviour, i.e. a smooth kinematic 
threshold of = Q^{x~^ — 1) = 4m|j, automatically. This method is cer- 
tainly appropriate provided is not too large. However, it does not sum 
potentially large logarithms in Q^/rrij^, and for this reason (and also because 
of the related lack of a heavy flavour parton distribution) is unsuitable for 

Recent 



at HERA 



there have been measurements of the charm structure function 
complementing the previous measurements by EMCn Also, the 
charm contribution to the total F2 at small x at HERA can be > 20%, and the 
NMC data, which contains only 5 — 10% charm, but has very small errors, is 
surprisingly sensitive to the treatment of charm. Hence, a global analysis of the 
most up to date structure function data must necessarily include a description 
of heavy flavours which is satisfactory over a wide range of x and Q^. 
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2 Theoretical Procedure 

We desire an approach which extrapolates smoothly from the FFNS at low 
to the ZM-VFNS at high Q^, maintaining the correct ordering in both 
schemes. In order to do this we first examine the relationship between the 
parton densities in these two schemes (putting /j,^ = for simplicity). One 
can show that the 4 flavour partons in the ZM-VFNS distribution are related 
to those in the 3 flavour FFNS by the equation 

ft{Q\^l) = AUQVml) ® fAQ'), (1) 

where the matrix elements Aka{Q'^ are calculable in perturbation theory, 
and contain logs in Q^/m^.^ Hence, the charm density is totally defined in 
terms of the light parton densities. Due to the formal equality of the two 
schemes for /m^ oo, in this limit we find a reciprocal relationship for the 
coefficient functions 

ClLiQVrnl) = ® AUQV^D- (2) 

So in this limit all the logs in ^^^(Q^/to^) can be factored into the matrix 
elements AkaiQ"^ /m'^), and absorbed into the definition of thejJi-flavour parton 
densities where they are summed automatically by evolutionu. 

For the scales < rn^ we assume the validity of the FFNS approach. For 
> we extend the approach outlined above to define a variable flavour 
number scheme (VFNS). We make precisely the same definition of the 4-flavour 
parton densities as above (thus ensuring that they evolve as massless partons), 
but in order to incorporate threshold effects correctly, i.e. guarantee exact 
equivalence to the FFNS, we have mass dependent coefficient functions. In 
particular we demand that for all > our coefficient functions satisfy 

CiJaiQV^c) = Cl^kiQVrnl) ® AkaiQVml). (3) 

Unfortunately, because there are more degrees of freedom in our coefficient 
functions than there are equations determining them, they are |iiot defined 
uniquely. For example at lowest order in as we have the equation □ 

<c?(QV-?) = Cl^"{QVml) ® A^,{QVml) + Cl^^'iQVml). (4) 

Since in the correct ordering of F2,c{x,Q^) only the zeroth order coefficient 
functions appear at LO, the first order coefficient functions appearing at NLO, 

"^At NLO (but not beyond) the charm distribution starts from zero at = and all 
other partons are the same in both schemes at this scale. 
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the way in which we choose to solve (Q), and higher order equations, affects 
our order-by-order expressions for _F2,c(a;, Q^). We remove this ambiguity by 
imposing the physical requirement that at the transition point = , where 
we go from a 3-flavour to 4-flavour description, dF2^c{x,Q^) / dh\Q^ is contin- 
uous. Imposing this to all orders (in the gluon sector) uniquely determines all 
coefficient functions and defines our VFNS.|^ Imposing this constraint, the LO 
coefficient function is given by the relatively simple expression 

5<f/(QV-^)/ainQ^ = cUfiQ^ml) ® P^^, (5) 

where P^g is the LO quark-gluon splitting function. This then automatically 

defines C2jg^{Q^ /ml) via (^. (^^^'^(Q^/jtIc) is more compficated, but numer- 
ically unimportant. The treatment for bottom quarks is essentially identical, 
and the procedure can easily be generalized to other processes. Further details 
can be found in El. 

Results 

Fig. 1^ shows the results of the NLO calculation of F2^c{x, Q^), illustrating how 
it smoothly extrapolates from the FFNS result at ml to the ZM-VFNS 

''We note that both our coefficierth functions and our definition of ordering is difl'erent to the 
alternative, currently LO, VFNSEl. 
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at high . We have also performed global fits using different prescriptions 
for heavy flavours and find that the best quality fit is indeed obtained using 
our VFNS approach Fig. || shows the prediction for i^2,c(a;, Q^) obtained 
from the parton distributions, resulting from our best fit to the total ^2(2;, Q^) 
(using TOc — 1.35GeV), for a variety of values of rric- Again one can see the 
smoothness and the good agreement with the full range of available data. 

Thus, our VFNS for heavy quark production is on a firm theoretical foot- 
ing, it is quite simple to implement, clearly has the behaviour we would expect, 
and matches well to data. Hence we encourage its use in analyses of structure 
functions. 



References 

f . Hf collaboration: C. Adloff et ai, Z. Phys. C 72 593 (1996). 

2. ZEUS collaboration: J. Breitweg et al, Phys. Lett. B 407 402 (1997); 
Paper N-645 presented at International Europhysics Conference on High 
Energy Physics, HEP97, Jerusalem 1997. 

3. J.J. Aubert et al, Nucl. Phys. B 213 31 (1983). 

4. M. Buza et al, Eur. Phys. J. C 1 301 (1998). 

5. M. Aivazis et al, Phys. Rev. D 5 3085 (1983). 



6. R.S. Thorne and R.G. Roberts, hep-ph/9709442, Phys. Rev. D 57 



(1998); R.S. Thorne and R.G. Roberts, Phys. Lett. B 421 303 (1998). 



4 



